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I present a calculation of structure functions at leading order which includes an 
unambiguous inclusion of the leading ln(l/a;) terms for each power of cis, and also 
the correct effects due to the mass of the charm and bottom quarks. I compare the 
results of fits to data to those obtained using conventional NLO in Os calculations, 
noting a clear preference for my approach, especially at small x. The predictions 
for both the charm structure function and Fl(x, Q^) using the two approaches are 
compared, the latter being much more discriminating. 



1 Introduction 

There has recently been a great deal of theoretical activity concerning the 
calculation of structure functions^ small x, driven by the vast increase in data 
in this region obtained at HERAEl'El. One of the main questions is whether one 
should include potentially important leading ln(l/2;) terms at high orders in a^, 
or simply order in powers of as alone. In a number of previous publications 
I showed that it is possible to include the leading ln(l/x) terms within the 
conventional renormalization group approach in an unambiguous manner at 
leading order El. The method of calculation is based on 3 points. 

1. The quantities one calculates to a given order are directly observable. Hence, 
I calculate in terms of the structure function at a particular scale (Q/) and its 
evolution away from this scale 

2. The leading-order expression for each independent part of a physical quan- 
tity begins at its lowest power of as, i.e. if a term ln"'(l/a;) first appears at 
order a" this is the leading order for this form of x dependence. 

3. The inputs for the structure functions are two flat nonperturbative functions 
for Fl and F2, convoluted with calculable perturbative contributions. These 
perturbative parts are determined by demanding that the expressions for the 
structure functions are invariant under changes of starting scale order by order 
in as, and are hence determined by the perturbatively controlled evolution. 

This leads to a leading-order-renormalization-scheme-consistent (LORSC) 
calculation of structure functions, which naturally combines the leading ln(l/ x) 
expansion with the more conventional expansion. When a global fit was 

"This results in expressions in terms of Catani's "physical anomalous dimensions'^. 



Table 1: Quality of fit using the LORSC{H) and the NLO-in-Os (MRST) approaches. 



x-range 


data points 


LORSC(H) 


MRST 


X > 0.1 


597 


688 


682 


0.1> X > 0.01 


385 


382 


377 


X < 0.01 


278 


219 


272 


total 


1260 


1289 


1332 



performed the quality using the LORSC calculation was superior to that using 
a conventional NLO-in-ag calculation, particularly at small x. 

2 Implementation of Heavy Quarks 

The main problem with this previous analysis was that it used a very naive 
prescription for heavy quarks, i.e. the charm and bottom quarks were both 
treated as infinitely massive below thresholds = mjj, and as massless ajbaye 
these thresholds. With the direct data on the charm structure function l3'l3i!l, 
and the fact that charm comprises ~ 20% of the total ^2(2:, Q^) at the lowest 
X values at HERA this is no longer sufficient. A method for including the 
heavy quark contributions in a manner which guarantees both smoothness 
at the threshold of = Q'^{x~^ — 1) = Arnjj and the correct supipation 
of ln{Q^/m'jj) terms was developed for the conventional approach llJ. The 
extension to the LORSC calculation is in principle quite simple, but in practice 
rather involved. Essentially it involves imposing matching conditions at = 
m'j^ such that the evolution is continuous at this point, but in terms of effective 
heavy quark coefficient functions and anomalous dimension^ which determine 
the heavy quark structure function and its evolution in terms of the light 
quark structure functions, rather than parton distributions. The details of 
this LORSC (H) calculation will be presented in a future publication. 

The quality of the LORSC (ED fit to a wide variety of structure function 
data (references can be found inEI) is compared to that for a NLO-in-ofg fit 
using an analogous treatment of charm, i.e. the recent MRST fitEI. The quality 
of the fit (using rric = 1.35GeV), in different x bins is shown in table |l|. The 
LORSC (H) is better overall than the MRST fit, being very slightly worse at 
high X, but considerably better at small x. 

We can also compare the results of the predictions from the two ap- 
proaches. Fig. shows the LORSC(H) calculation of F2^cix,Q^) (using 
nic — 1.35GeV and rric = 1.5GeV), and also the prediction from the MRST fit 
(where rric = 1.35GeV). Although in principle the two approaches could give 

''These involve the leading ln(l/a;) heavy quark coefiicient functions^. 
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' Q'(GeV-) "' 

Figure 1: Comparison of F2,c{x, Q^) with data using the LORSC(H) fit with rria = f .35GeV 
and f.SGeV and the MRST fit {rUc = 1.35GeV). 

different predictions for F2,c{x,Q'^), in practice they are rather similar when 
the same value of nic is used. However, while the MRST fit is very sensitive to 
rric, becoming worse quickly when it increases above 1.35GeV, the LORSC(H) 
fit is almost unchanged in going from 1.35GeV to l.SGeV, and provides much 
more freedom in F2^c{x,Q^). A larger difference between the approaches is 
observed when comparing the predictions for Fl{x, Q^), as seen in fig. ||. We 
note that both predictions for Fl(x, Q^) are significantly lower than when us- 
ing the previous treatment for charm, since in the correct approach the heavy 
quark contribution to Fl(x, Q^) is strongly suppressed until 3> rn^. 

3 Conclusion 

It appears as though the inclusion of leading ln(l/x) terms in a consistent man- 
ner significantly improve the comparison to structure function data at small 
X, highlighting the shortcomings of a NLO-in-a^ calculation in this region. 
The calculation of the kernel _£Dr_the NLO-in~ln(l/x) gluon Green's function 
has recently been completed llilllj, and is very suggestive that the NLO-in- 
ln(l/a;) corrections are extremely large. Nevertheless, further understanding 
seems necessary before it is known precisely how this new result relates to a 
consistent calculation of structure functions. It may indeed be true that a RSC 
expansion scheme is not really convergent. However, it is already known that 
this is true for the conventional as expansion: the predictions for Fl{x,Q^) 
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b igure 2: Prediction for Fl(x,Q^) in the HERA range for the LORSC(H) fit (lower curves) 
and for the MRST fit (upper curves). 



being hugely different at small x at LO and NLO. Hence, I propose that the 
sucess of the inclusion of leading ln(l/a;) terms in the correct manner is telling 
us something important about the true physics at small x. 
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